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^ ; ABSTRACT 

Thanks to Hubble and Chandra telescopes, some of the large scale jets in ex- 
tragalactic radio sources are now being observed at optical and X-ray frequencies. 
For the FR I objects the synchrotron nature of this emission is surely established, 
although a lot of uncertainties - connected for example with the particle accel- 
eration processes involved - remain. In this paper we study production of high 
energy 7-rays in FR I kiloparsec-scale jets by inverse-Compton emission of the 
synchrotron-emitting electrons. We consider different origin of seed photons con- 
tributing to the inverse-Compton scattering, including nuclear jet radiation as 
well as ambient, stellar and circumstellar emission of the host galaxies. We dis- 
cuss how future detections or non-detections of the evaluated 7-ray fluxes can 
provide constraints on the unknown large scale jet parameters, i.e. the magnetic 
field intensity and the jet Doppler factor. For the nearby sources Centaurus A 
and M 87, we find measurable fluxes of TeV photons resulting from synchrotron 
self-Compton process and from comptonisation of the galactic photon fields, re- 
spectively. In the case of Centaurus A, we also find a relatively strong emission 
component due to comptonisation of the nuclear blazar photons, which could be 
easily observed by GLAST at energy ~ 10 GeV, providing important test for the 
unification of FR I sources with BL Lac objects. 
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Subject headings: radiation mechanisms: nonthermal — galaxies: jets — galaxies: 
individual (Centaurus A, M 87) — gamma rays: theory 

1. Introduction 

7-rays are detected from many galactic and extragalactic sources, including pulsars, 
supernova remnants, gamma ray bursts and blazars. It is believed that multiwavelength 
emission of the latest type of the mentioned objects is due to small scale highly relativistic 
jets, which are formed near active galactic nuclei of host galaxies and are inclined at small 
angles to the line of sight (e.g., Bregman 1990). High observed bolometric luminosities of 
blazars (typically 10 44 — 10 48 erg/s) are often dominated by 7-ray emission, what enabled 
to observe many of them by EGRET in e 1 ~ 0.1 — 10 GeV energy range (von Montigny et 
al. 1995; Hartman et al. 1999). Some nearby low-luminosity blazars, like the best known 
Mrk 421 and Mrk 501, are also confirmed to be sources of very high energy (VHE) 7-rays 
with e 7 > 100 GeV (Punch et al. 1992; Quinn et al. 1996, respectively). Detection of such 
VHE radiation by ground-based instruments was possible due to development of imaging 
atmospheric Cherenkov telescopes, IACTs (see, e.g., Aharonian 1998; Volk 2003). 

The observed broad-band blazar emission allows one to determine physical parameters 
of nuclear AGN jets. 7-ray observations are crucial in this respect, providing important 
constraints on the involved particle acceleration processes, or on the characteristic spatial 
scales for the blazar phenomenon (e.g., Sikora and Madejski 2002; Kino et al. 2002). However, 
observations at VHE range are complicated because propagation effects are important for the 
high energy photons created at cosmo logical distances (Nikishov 1962; Gould and Schreder 
1966; Stecker et al. 1992), and also, obviously, because of a low photon statistics in 7-ray 
telescopes. Among the other issues, one should mention here a case of VHE emission from 
Mrk 501, which extends up to energies e 7 > 10 TeV, which seems to be in conflict with some 
models of cosmic infrared background (CIB) radiation (see, e.g., Protheroe and Meyer 2000; 
Aharonian 2001). Detailed studies of cosmic background photon fields and their interactions 
with 7-rays during their propagation, as well as future 7-ray missions, will possibly answer 
some of the questions connected with production of VHE radiation in blazars. 

Contrary to blazar parsec-scale jets, their large scale counterparts extending from a few 
to a few hundreds of kiloparsecs from active galactic nuclei were usually studied only at 
radio frequencies (Bridle and Perley 1984). Recently, however, Hubble (HST) and Chandra 
telescopes gathered new information about optical and X-ray emission of some of these 
objects. For several hundreds of known radio jets, only about 20 are observed at optical 
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frequencies 1 . Most of them are relatively short and faint, with only a few exceptions like the 
one in M 87 allowing for detailed spectroscopic and morphological studies (Meisenheimer 
et al. 1996; Perlman et al. 1999, 2001). Surprisingly, the large scale jets can be also very 
prominent in X-rays. Up to now, more than 20 jets were detected by Chandra within the 
1 — 10 keV energy range 2 , although the nature of this emission is still under debate (e.g., 
Tavecchio et al. 2000; Bai and Lee 2001; Celotti et al. 2001; Aharonian 2002; Carilli et al. 
2002; Harris and Krawczynski 2002; Stawarz and Ostrowski 2002; Sahayanathan et al. 2003; 
Tavecchio et al. 2003). In general, both optical and X-ray observations question in many 
aspects the standard models for the extragalactic jets emission. 

As the large scale jets in radio galaxies are confirmed sources of optical and X-ray 
photons, one can ask if they can be also observed at 7-ray energies. Such observations could 
impose significant constraints on jet parameters (in analogy to the blazar sources), as well 
as for the unification scheme of radio-loud AGNs. In the present paper we make predictions 
regarding production of high energy 7-rays in FR I jets by the inverse-Compton scattering 
in a framework of one electron population model. In discussion, we refer to the objects 
observed by Chandra, as the X-ray observations of these recently widely studied sources 
impose important constraints on the jet synchrotron emission, crucial in determining the 
high energy photon flux. We find, that for the jet parameters inferred from radio-to- X-ray 
observations, only nearby Centaurus A and M 87 jets are expected to be observed at GeV 
- TeV energies by near future 7-ray missions (cf. discussion in section 3). However, it is 
possible that more distant objects will be also detected, if the unconstrained jet parameters 
(like magnetic field and kinematic factors) differ from typically assumed values (sections 
2.3.1 - 2.3.3). 

The already mentioned propagation effects can seriously affect the observed VHE spectra 
of distant objects, as interaction of VHE 7-rays with the CIB photon field is expected to 
result in their efficient absorption due to photon-photon pair creation. CIB radiation is 
created by infrared emission of forming stars in early cosmological epochs and the following 
reprocession by ambient dust. The origin, evolution and detailed energy distribution of CIB 
are still uncertain, mostly because of observational difficulties (for a recent review see Hauser 
and Dwek 2001). The existing controversies, particularly important in a TeV-absorption 
range (at u C ib ~ 10 13 Hz) are however not crucial for our discussion, as we limit ourselves 
to relatively close object. For illustration, we refer to work by Kneiske et al. (2002, and 
references therein) who, among the others, studied opacity of the Universe to 7-rays and 



1 see homepage http://home.fnal.gov/~jester/optjets by S. Jester and references therein 
2 see homepage http://hea-www.harvard.edu/XJETS by D. Harris and references therein 
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evolution of the diffusive infrared-to-ultraviolet background radiation field. 

Below, in section 2, based on optical and X-ray observations we reconstruct and comment 
'typical' electron energy distribution of FR I Chandra jets. We estimate different photon 
fields contributing to the inverse- Compt on scattering of these electrons, and next evaluate 
and discuss the resulting fluxes and photon break energies in the Thomson regime. Section 
3 contains detailed analysis of the 7-ray emission from the nearby sources Centaurus A and 
M 87. Final conclusions are presented in section 4. 

2. Emission from the large scale jets in FR I sources 

2.1. Synchrotron emission: radio-to- X-rays 

Chandra X-ray observatory detected about 10 jets displaying the FR I large scale mor- 
phology. They include the ones in radio galaxies M 87 (Marshall et al. 2002; Wilson and Yang 
2002), Centaurus A (Kraft et al. 2000, 2002), 3C 129 (Harris et al. 2002a), 3C 31 (Hardcastle 
et al. 2002), PKS 0521-365 (Birkinshaw et al. 2002), 3C 270 (Chiaberge et al. 2003), M 84 
(Harris et al. 2002b), 3C 66B (Hardcastle et al. 2001), B2 0206+35 and B2 0755+37 (Worrall 
et al. 2001). All of these objects are located relatively nearby, with a distance ranging from 
3.4 Mpc (Centaurus A) to more than 300 Mpc (PKS 0521-315 at redshift z = 0.055). All of 
them share many spectral and morphological similarities briefly summarised below. 

X-ray jets in FR I radio sources are usually quite short (projected length ~ 1 — 4 kpc). 
Typically, they are composed of diffusive knots with a spatial scale it! ~ 0.1 kpc, although 
a strong inter-knot emission is also sometimes present. X-ray jet morphology corresponds 
roughly to the radio morphology, and - in cases of M 87, PKS 0521-315, 3C 66B, 3C 31 
and B2 0755+37 observed by HST - to the optical one. A noted difference between these 
pictures are spatial offsets of some knot maxima as measured at X-rays and at radio/optical 
frequencies (up to ~ 0.008 kpc in M 87, ~ 0.08 kpc in Centaurus A and ~ 0.2 kpc in 3C 
66B). Also, X-ray jets (knots) seem to be narrower than their radio/optical counterparts (M 
87, Centaurus A). Except of the weakest and the smallest objects (3C 270 and M 84), the 
observed X-ray luminosities of the discussed jets are L x ~ 10 39 — 10 42 erg/s. 

In most cases, multiwavelength observations of knot regions in FR I Chandra jets allow 
one to construct broad band radio-to-X-ray spectral energy distribution. All of the inferred 
knot spectra seem to be similar, and in addition difficult to be explained in a framework of 
standard jet models. A spectral index of the radio emission (defined in a way S u oc v~ a , 
where S v is the energy flux spectral density) is always close to or ~ 0.6. The optical emission 
obviously belongs to the synchrotron continuum, with the spectral index steepening between 
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infrared and optical frequencies at about Ao ~ 0.7 (3C 66B, PKS 0521-315, M 87). This 
agrees well with the X-ray spectral index measured by Chandra, which, for FR I jets, is on 
average ax ~ 1-3 (ranging from ~ 1.1 in 3C 31 to > 1.5 in some knots of M 87 and Centaurus 
A). Thus, one can generalize, that in the discussed objects the radio-to-optical power-law 
slope is always less than unity (ctRo ~ 0.5 — 0.8), while the optical-to- X-ray power-law slope 
is always much larger (aox ~ 11 — 1-5). Such spectral behaviour seems to be universal 
for all FR I X-ray jets, with a spectral break (br) - where most of the synchrotron (syn) 
power is emitted - placed usually around v S yn,br ~ 10 14 Hz or, eventually (M 87), at higher 
frequencies. 

The steep X-ray spectral indices of FR I large scale jets exclude an inverse-Compton 
scattering as a mechanism responsible for the observed X-ray emission. Instead, radio-to- 
X-ray spectral properties indicate synchrotron origin of detected 1 — 10 keV photons. If 
this is the case, than in the discussed objects unknown acceleration process operates, which 
results in characteristic and universal for FR I jets electron energy distribution which can 
be approximated as a broken power-law: 



where 7 is the electron Lorentz factor 3 , Aa is the spectral break in the synchrotron continuum 
and p = 2 an + 1 ~ 2.2. Optical and X-ray observations play a crucial role in determining 
parameters of the above electron spectrum, i.e. its normalisation (for the assumed energy 
equipartition between electrons and the jet magnetic field) and the actual values of % r 
and ■jmax, which in turn determine the inverse-Compton 7-ray jet radiative output. For non- 
relativistic jet velocities, a typical v syn ^ r ~ 10 14 — 10 15 Hz and an equipartition magnetic field 
B ~ 10 -4 G, one gets 7 br ~ 10 5 — 10 6 and 7 ma;E > 10 7 . Discussion of the acceleration processes 
creating such high energy electrons in the FR I jets is beyond the scope of this paper. Let 
us note however, that it is possible for 71^(7) to posses a more complicated form than the 
one given by equation 1, involving for example spectral pile-ups at the maximum electron 
energies due to efficient electron acceleration and rapid radiative cooling (Dermer and Atoyan 
2002; Stawarz and Ostrowski 2002). In either the observed X-ray luminosity of the 

FR I kpc-scale jets is lower than the optical luminosity, most of the electrons energy is 
concentrated in the range near the electron break Lorentz factor % r . 



3 Hereafter we follow the notation with the primed quantities measured in the jet comoving frame and 
the bare ones if given in the observer rest frame (neglecting cosmological corrections). However, we do not 
prime the magnetic field induction B as well as electron Lorentz factors 7, noting instead that they always 
refer to the emitting plasma rest frame. 




for 70 < 7 < Ibr 

for 7fe r < 7 < 'Jmax 



(1) 



- 6- 



It is interesting to compare v sy n,br characteristic for the kpc-scale FR I jets with critical 
(cr) synchrotron frequency v syncT resulting from the interplay between electron synchrotron 
cooling and dynamical evolution of the emitting region. The time scale for the former 
process can be estimated as t' loss ~ t/ItIsj/td where mc 2 \ < y\ S yn — § C(T t^b7 2 is a mean rate 
of electron energy losses due to synchrotron emission, and U' B = B 2 /87i denotes comoving 
jet magnetic field energy density. The second time scale is simply t' dyn = td yn /F ~ r/cT, 
where T = (1 — /3 2 )™ 1 / 2 is the jet bulk Lorentz factor, (3c is a jet bulk velocity and r is 
the distance of the kpc-scale structure from the jet base at r C r from the active galactic 
nucleus. A critical electron Lorentz factor is defined by an equality t' dyn ~ t' loss {j = 7 cr ). 
Thus, assuming initial injection of the power-law electron energy distribution at r , one 
expects that at the distance r the electron spectrum steepens for 7 > 7 cr because of the 
efficient radiative cooling of these electrons. The observed critical synchrotron frequency is 
then v ayn ^ cr = v' syncr b~ oc "f 2 cr B5, where 

6 = r(l-/?cos0) (2) 

is the jet Doppler factor and 9 is the jet inclination angle to the line of sight. Assuming 
that the jet magnetic field scales with the distance as B = B (r /r), what conserves the 
Poynting flux L' B oc R 2 U' B in the expanding jet with an opening angle f ~ const and a 
radius R ~ ipr, one obtains 

v cr ~4-10 14 rx5r 2 Hz , (3) 

where we put r = 1 pc, B — 0.1 G and r 1 = r/1 kpc. Note, that with the above values of 
Bq and ro characteristic for the blazar sources, the kpc-scale jet magnetic field is expected 
to be B ~ 10~ 4 G, consistently with the equipartition value. Hence, in a case of a non- 
relativistic kpc-scale jet velocity, v syn ^ cr ~ v S yn,br ~ 10 14 Hz. However, the spectral break 
resulting from the considered process is Ao; cr = 0.5 (Kardashev 1962), i.e. slightly less than 
the one required for the FR I Chandra jets. Possibly, the observed steep spectral break is a 
signature of the electron synchrotron cooling in the spatially inhomogeneous magnetic field 
(Cavallo et al. 1980; Coleman and Bicknell 1988). Note also, that in the cases of relativistic 
jet bulk velocities, the observed value of v syn ^ CT can be significantly higher than 10 14 Hz (cf. 
the case of M 87, section 3.2). 

Radio observations of FR I jets indicate moderate or even weak beaming (e.g., Laing et 
al. 1999). However, one cannot exclude possibility that the bulk Lorentz factors of kpc-scale 
jets in weak radio galaxies are of order of a few. This idea is supported by EST observations 
of superluminal motions and Chandra detection of significant X-ray variability at kpc-scales 
in M87 (Biretta et al. 1999; Harris et al. 2003, respectively). In fact, a small number of 
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detected optical and X-ray jets in radio galaxies as compared to the number of known radio 
jets suggests that relativistic beaming effects can play an important role in these sources 
(Sparks et al. f995; Scarpa and Urry 2002; Jester 2003). Including a relativistic correction, 
the electron break Lorentz factor is 

*~ 5 ' 105 fey) • < 4 > 

where v 8y n,u = v syn ^ r j\^ Hz and _B_ 4 = 5/10 -4 G. 



2.2. Radiation fields within the kpc-scale jet 

Let us consider a kpc-scale (r ~ 1 — 4 kpc) relativistic jet with the bulk Lorentz factor 
T, inclined at an angle 9 to the line of sight. The energy density of the jet magnetic field, as 
measured in the emitting region rest frame, is 

B 2 

U' B = 4 • lO" 10 ^ erg cm" 3 . (5) 

8 7T 

As mentioned above, the equipartition value inferred from the synchrotron emission of FR 
I knot regions is typically B eq ~ 10~ 4 G, and it corresponds to the observed synchrotron 
jet luminosities L syn ~ 10 40 — 10 42 erg/s. The observed luminosity is related to the total 
emitted synchrotron power L' syn by the relation 

L S yn = g c j/ m b(T, 9) L' syn , (6) 

where 

a (T9)-S 63/T M (7) 

9a/um-[ 5 4 {mb) (7) 

in cases of a continuous jet (cj) or a moving single radiating blob (mb), respectively (Sikora 
et al. 1997, see also Appendix A) 4 . Below we analyse these both possibilities as limiting 
models for the kpc-scale knots. We also consider synchrotron luminosity at a given break 
frequency v syn ^ r , hereafter denoted as [vL u ] syn ^ r , rather than the bolometric one L syn , with 
transformations between \yLj\ syn ^ r and \y' ' L'^syn^ the same as given in equations 6-7. 
For the broken power-law synchrotron spectrum assumed here, with a ~ 0.6 and Aa ~ 0.7 
(cf. equation 1) one has a bolometric correction L syn /[uL u ] syn ^ br 6. 



4 Notc, that as discussed in Appendix A, in both cases the equipartition magnetic field measured in the 
emitting region rest frame is related to the equipartition value computed for no beaming by the relation 

B eq = B eq ^ $=1 S~ 5 / 7 . 
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In addition to the synchrotron emission extending to the X-rays, the FR I jets produce 
also the high energy radiation by inverse- Compton scattering of ambient photon fields, like, 
for example, the cosmic microwave background (CMB) radiation or the synchrotron emission 
produced by the jet itself. At the kiloparsec distances from the galactic center another 
important sources of ambient photons are the active galactic nucleus and the host galaxy. 
The AGN radiation, produced at small distances from the central region, illuminates the 
large scale jet almost exactly from behind. This radiation consists of isotropically emitted 
component connected with thermal gas and/or dust heated by the central source (the narrow 
line region and the dusty nuclear torus) plus an anisotropic blazar-like emission due to small 
scale highly relativistic nuclear jet. The narrow line emission and the radiation of dusty 
nuclear tours, typically weak or even absent in FR I radio galaxies (e.g., Chiaberge et al. 
1999), are always negligible at distances r > 1 kpc from the active nucleus (section 2.2.2, 
see also Celotti et al. 2001). Instead, as discussed below, comptonisation of starlight and 
extended dust emission of the host galaxy and, in some cases, comptonisation of the 'hidden' 
blazar radiation can dominate 7-ray output of the FR I large scale jets. 



2.2.1. Synchrotron photons 

For a simple evaluation of the synchrotron photons energy density in the jet comoving 
frame, we consider an approximately cylindrical knot region, with a radius R < 0.1 kpc and 
a deprojected observed length I ~ R. In the case of a continuous jet one has V = IT > R, 
and therefore 

2rcRl'c ~ 2vRlc P ™ ' 10 R\ P ' (§) 

where [vL u ] synt ±2 = [vL v } S y n ^ r /lQ A2 ergs -1 and = i?/0.1kpc. If the emitting kpc-scale 
jet region is in fact a moving source, then I' = I / 5 = R / 5, and the jet comoving energy 
density of the synchrotron emission is 

The jet comoving synchrotron break frequency is v' syn hr = v S yn,brl^- Note, that for the 
approximately constant ip ~ 0.1, the jet radius at a distance r ~ 1 kpc from its base is 
expected to be R ~ 0.1 kpc, consistently with the value considered here. 
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2.2.2. Radiation of the active nucleus 

At first, let us consider the blazar-like (bl) emission of FR I active nuclei, which - 
accordingly to the unification scheme (e.g., Urry and Padovani 1995) - belong to the low- 
luminous (BL Lac) blazar subclass. Due to relativistic effects, this radiation is strongly 
beamed into the cone with an opening angle 9 bl ~ where T bl is the nuclear jet bulk 

Lorentz factor. In order to evaluate energy density of this emission in the kpc-scale jet rest 
frame, we assume that it 'enters' into the considered emitting region directly from the jet 
base. In other words, we neglect any possible misalignment between the pc-scale and the 
kpc-scale jets. The respective energy density of the blazar photons is then 

^~4^(2Tp ' (10) 

where L b i(0) is the isotropic luminosity of the blazar emission evaluated by the stationary 
observer located at 9 = (Appendix B). With the spatial scales of the kpc-scale jet emitting 
region much larger than the characteristic nuclear scales (R, r ^> r ), the discussed source 
of the seed photons for the inverse- Compton scattering should be considered as a stationary 
one. A typical short time scale of blazar variability justifies this statement. Therefore, 

L " (0) = Lu -ig^r Lu { 4tJ ■ 

where 5 b i = [T b i (1 — f3 b i cos#)] -1 and L b \ is the isotropic luminosity of the considered sta- 
tionary blazar source measured by the observer located at the angle 6 to the jet axis. Hence, 

Ubl 2^1p 5 - 6 ' 10 ^PTp ergcm ' (12) 

where we put L W 42 = L w /10 42 erg s" 1 (cf. equation 1 in Celotti et al. 2001). 

Typically, for the high energy peaked BL Lacs (HBLs) the observed synchrotron lumi- 
nosity is 10 44 — 10 46 ergs _1 , while for the low energy peaked BL Lacs (LBLs) it is in the 
range 10 45 — 10 47 ergs _1 (Fossati et al. 1998). An another subclass of blazar sources - the 
flat spectrum radio quasars (FSRQs) - have even higher observed luminosities. However, in 
a framework of the unification scheme FSRQs correspond to the FR II sources. The critical 
observed synchrotron frequency of blazar emission, which anticorrelates with the blazar syn- 
chrotron luminosity, is usually 10 15 — 10 17 Hz for HBLs and 10 13 — 10 15 Hz for LBLs (again 
Fossati et al. 1998). Below, we take L b \ ~ [uL u ] b i tbr , where the blazar synchrotron break fre- 
quency in the jet comoving frame is v' bl br ~ i^>(0)/r, and v b i, br (0) = v UM (2T b i/5 b i). We 
also assume, that the blazar source illuminates uniformly the whole kpc-scale jet emitting 
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region. Note, that for the typical T b i ~ 10, the half-angle of the blazar emission cone is equal 
to the assumed jet opening angle, 6u ~ ip ~ 0.1. 

As mentioned previously, the isotropic component of the nuclear (nucl) emission, i.e. 
radiation of narrow line region or dusty nuclear torus, are weak or even absent in FR I 
radio galaxies (isotropic luminosity L nud < 10 42 erg/s). For the observer located at the 
distance r from the active center and moving with the Lorentz factor T along the jet axis, 
energy density of the considered circumnuclear isotropic emission scales in the same manner 
as the energy density of the blazar radiation - in both cases the appropriate relations are 
L n ud/^r 2 c(2T) 2 and the one given in equation 10, respectively. Therefore, as 
Lnud <^ L b i(0), the isotropicly emitted nuclear photon fields are negligible as compared to 
the blazar emission. 



2.2.3. CMB and galactic photon fields 

AGNs are hosted by elliptical galaxies, whose extended optical and infrared emission 
contribute significantly to the total radiation intensity at the kiloparsec distances from the 
galactic centers. 

Ellipticals are known to contain ~ 1 M & stars evolving from the main sequence through 
the red-giant phase, hot interstellar gas heated to temperatures ~ 10 7 K and radiating at 
X-rays, and, finally, dust prominent at infrared frequencies. Photospheric emission from 
cool giant stars, which are concentrated within a few hundreds of parsecs from the galactic 
center (i.e. within a 'central bulge' or a 'galactic core'), dominates bolometric luminosities 
of the discussed galaxies and peaks in the range ~ 1 — 3/im (e.g., Knapp et al. 1992). Total 
bolometric luminosities of giant ellipticals can be as high as ~ 3 • 10 11 L Q ~ 10 45 erg/s. 
Emission at wavelengths 5 — 20 is analogously distributed as the starlight and is most 
likely connected with dusty winds accompanying the red-giant stars and forming their dusty 
circumstellar envelopes (Knapp et al. 1992; Tsai and Mathews 1995). Radio galaxies are in 
general more luminous at infrared wavelengths as compared to the normal ellipticals (by a 
factor of 2 - 3 relative to the starlight emission), with ~ 10 /im luminosity being on average 
two orders of magnitude lower than the bolometric one (Knapp et al. 1990). In addition to 
the circumstellar dust emission, ellipticals exhibit also an excess at far-infrared frequencies, 
peaking around 60 — 100 /i m. If the cold dust responsible for production of these photons is 
optically thin, it must occupy regions that extend far from the galactic core (a few times the 
core radius), in order to avoid overheating by the starlight (Tsai and Mathews 1996). An 
alternative explanation can be provided by concentration of the cold dust in disks and lanes 
self-shielded from the starlight and located closer to the galactic center. The far-infrared 
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excess is especially prominent in radio galaxies (Golombek et al. 1988; Knapp et al. 1990). 
However, as the origin and spatial distribution of this radiation is unknown, and also as the 
far-infrared emission shows significant variations from galaxy to galaxy, we do not discuss 
its contribution to the kpc-scale jet inverse-Compton emission. 

In Appendix C, following Tsai and Mathews (1995), we estimate the energy density 
of the starlight (star) photons to be U star ~ 10 -9 erg/cm 3 at the distance ~ 1 kpc from 
the center of the typical giant elliptical galaxy. Assuming for simplicity an approximately 
isotropic distribution of the stellar emission at this scale, in the jet comoving frame one has 

U'star ~U star T 2 ~\^ % T 2 erg cm" 3 (13) 

(Appendix B), with a characteristic starlight frequency v' star ~ 10 14 T Hz. The galactic 
dust emission (dust) is distributed analogously to the stellar radiation and hence its energy 
density in the jet rest frame is also amplified approximately by a factor T 2 . In this paper we 
take 

U' dust ~ 0.01 • U' star ~ 10- 11 T 2 erg cm" 3 (14) 

at v' dust ~ 3 • 10 13 T Hz. Finally, for redshifts z € 1 an analogous energy density of the 
blackbody CMB is equal to 

U' CMB = aT^ MB T 2 ^A-10- 13 T 2 erg cm" 3 , (15) 

where a = 7.53 ■ 10~ 15 cgs and the observed CMB temperature is Tqmb — 2.7 K. The 
appropriate characteristic CMB photon frequency is v' CMB ~ 2 ■ 10 11 T Hz. 



2.3. Inverse-Compton emission: 7-rays 

With the evaluated photon fields in the emitting region rest frame, U' seed , one can esti- 
mate the observed break luminosity of the appropriate inverse-Compton (ic(seed)) emission 
in the Thomson regime as 

\vLi\ic(seed)M ~ f±/iso(T, 9) —— [uL^ syn ^ b r (16) 

Presence of an additional factor f±/i SO (^, 6) is connected with the possible anisotropy of 
the external radiation fields within the jet comoving frame, in cases of its relativistic bulk 
velocities and/or non-isotropic photon distribution in the galaxy rest frame. As shown in 
Appendix D, 

f (5/V) 2 [(l + ^)/(l + P)} 2 (+) 
f±iUr.9)={ \(5T) 2 [(l-^)(l + (3)f (-) (17) 

1 (iso) 
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respectively for the cases when the seed photons are distributed isotropically in the galaxy 
rest frame near the kpc-scale jet (+), when they illuminate the jet exactly from behind (— ), 
or when they are isotropic in the jet comoving frame (iso). Luminosity [^^(seed),^ is 
related to the observed flux of the inverse- Compton emission, [vS v ] ic ( seed ^br-> by the relation 

Lrrl _ \yi-<v\ic{seed)M / 1fi \ 
[VOv\ic(seed),br ~ ^ ^ 2 , 

where D is a distance to the source. The appropriate observed inverse- Compton break 
frequency due to 7^ electrons scattering the seed photons with the observed characteristic 
(break) frequency u seed is 

^ic(seed),br 

h±/iso(T, 6) 

T&r V seed i 

(19) 

where 

r 5 2 (i + ^)/(i + /3) (+) 

h±,iso(r,0)=\ 5 2 (l-fi)(l + P) (-) (20) 
[ 4/3 {iso) 

(Appendix E). Assuming that the scattering proceeds in the Thomson regime, 

lbr<eed<— , (21) 

the inverse-Compton spectrum for v ic ( seed ) < u ic ^ seed ^ br is a power-law with a spectral index 
corresponding to the radio-to-optical synchrotron continuum, a 7 ~ {p — l)/2 (equation 1), 
while for v^eed) > ^ic(seed),br it is expected to steepen by Act. The high energy cut-off in 
the observed 7-ray emission can result from the cut-off in the electron energy distribution 
(imax), from entering the Klein-Nishina (KN) regime or due to external absorption on the 
CIB radiation. 

Below, we evaluate fluxes and break energies of the inverse-Compton emission on photon 
fields described in the previous section. At first we discuss the synchrotron self-Compton 
radiation, i.e. the inverse-Compton scattering of the synchrotron photons produced by the 
same electron population (ic = ssc; seed = syn). Next, we analyse the external-Compton 
{ic = ec) process starting from comptonisation of the 'hidden' blazar emission [seed = bl). 
Then we discuss the scattering of CMB photons and different radiation fields of the host 
galaxy {seed = star, dust, CMB). 
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2.3.1. Synchrotron self-Compton emission 

From the equations 8-9 and 16 - 18 with /j SO (r, 9) = 1, one obtains the observed SSC 
energy flux for the two considered cases of a continuous jet and a moving blob 

-11 W L v] 2 syn,42 f S~ 3 (cj) erg 



KUtr - 1-2 • io- n ^ 2 T'" ;_ 4 ; 7, (22) 

Bi 4 R± 1 D( { d 4 (mo) scm 2 

respectively, where Di = D/10 Mpc. Equations 19 - 20 with /ij SO (r, 6 1 ) = 1 give the observed 
SSC break energy 

£ssc> ~ 1.5 • 10 11 ^ 1 eV . (23) 
The Thomson regime condition for the discussed process reads as 

v S yn,u<2B 1 _ / ?5 . (24) 

The SSC emission is determined by the observed parameters [vL u ] synthr and v syn ,br, plus 
the emitting region linear size R, the (unknown) magnetic field B and the kinematic factor 8. 
Typically, for the FR I Chandra jets one has [vL u ] syn ^ 2 < 1, R-i ~ 1 and the equipartition 
value (as computed for nonrelativistic bulk velocities) 5_4 ~ 1. Hence, for v syn ^ ~ 1 and 
the usually considered 8 ~ 1, one can put an upper limit {vS v } SSCt b r < 1CT 11 -Df 2 erg/s cm 2 at 
the critical photon energy e ssc ^ r ~ 0.1 TeV. Future IACT systems will be able to eventually 
detect such emission 5 from the sources at the distances D < 100 Mpc. Higher SSC flux than 
discussed above can be produced for a given [h / L l/ \ syni i )r in the case of a strong departure from 
the equipartition, B <C B eq , and/or for a de-beaming of the jet emission, <Jc 1. However, 
for decreasing B and 8 the effects connected with the KN regime are expected to become 
important (unless v S yn,u <C 1, equation 24), decreasing the efficiency of the SSC process. 



2.3.2. Comptonisation of the hidden blazar emission 
From the equations 12 and 16 - 20, one obtains the observed EC(bl) energy flux 

r„o i in m-i3 \ vL v\biAi [vL v \ synA2 T\i f- erg 
[^W0^~l-2-10 B^rlD^l P W ' (25) 



5 with the detector sensitivity > 10 13 erg/scm 2 for 5a detection threshold and 100 h observation (Aha- 
ronian 1998) 
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and the observed EC(bl) break energy 

r> r> 1 nil V bl, 14 V syn 14 Ttf Jl- . . 

e ec (bi),br ~ 2.2- 10 11 ' y : F eV . (26) 

The Thomson regime condition for the considered process reads as 

r H v bl , u < r 5^ b x H v; y l( 2 u . (27) 

The EC(bl) emission is determined by the jet parameters [^£i/] SJ/n ,&r and v syn ,br plus the 
distance from the blazar source r, the unknown magnetic field B, kinematic factors T and 9 
and, finally, by properties of the blazar core [vL^i^, fu,br and T^. For the usually discussed 
8 ~ 1 and 9 > 45°, typical v syn ^ ~ 1 and the equipartition value _B_ 4 ~ 1, one expects the 
observed EC(bl) break energy to be £ ec (HBL),br ~ 1 — 100 TeV in the case of comptonisation of 
HBL radiation (bl = HBL), and £ ec (LBL),br ~ 0.01 — 1 TeV in the case of LBL-like core emis- 
sion (bl = LBL). As the synchrotron break frequency of the blazar emission is anticorrelated 
with the blazar synchrotron luminosity (see section 2.2.2), one can put upper limits on the 
observed EC(bl) fluxes {vS u ] ec (HBL),br < 10~ n _D]~ 2 erg/scm 2 for £ e c(HBL),br ~ 1 TeV (taking 
a rough estimate L HB L,max(0) < 10 46 erg/s), and [vS v \ ec ( LB L),br < 10~ 10 £>f 2 erg/scm 2 for 
£ ec (LBL),br ~ 10 GeV (with LL B L,max(0) < 10 47 erg/s). In addition, the EC(HBL) emission 
is expected to be significantly decreased due to the KN effects, contrary to the EC (LBL) 
radiation. For the discussed standard jet parameters the latter one could be observed by 
GLAST® from sources with D < 100 Mpc. Similarly to the case of the SSC process, for a 
given [vL v ] syn ^ r the subequipartition magnetic field and/or large jet inclination increase the 
expected value of [^S^ec^o,^- 

2.3.3. Comptonisation of the galactic photon fields and CMB radiation 
From the equations 13 and 16 - 20, one obtains the observed EC(star) energy flux 

[„(J1 oi 1 n-10 [vL^syn, 42 f p 2 er g / 0Q \ 

[^^u\ec(star),br ~ 2.1 • 10 r2 n2 7+1 — , (28) 

and the observed EC (star) break energy 

e«(^r), 6 r~l.l-10 11 ^^ eV . (29) 

-D-4 



6 accordingly to the 10 GeV GLAST sensitivity ~ 10 12 erg/s cm 2 for 5cr detection threshold in 1 year 
all-sky survey (Aharonian 1998) 
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The considered process proceeds within the Thomson regime for 



Vsyn, 14 < 6 5_ 4 — . (30) 



The EC(star) radiation is determined by the jet parameters [vL v ] synthr , v syn ^ hr , B and 
5, as well as by properties of the host galaxy (L star , etc). For non-relativistic jets, typical 
v S yn, 14 ~ 1, [^i/] syn, 42 < 1 an d -B-4 ~ 1, the break energy of the considered emission is 
£ec(star),br ~ 0.1 TeV, and the observed energy flux is [vS v \ ec ( sta r),br < 10~ 10 erg/scm 2 . 
However, this relatively strong radiation can be significantly decreased in the case of rela- 
tivists jet velocities and large jet inclinations to the line of sight due to both the KN and 
the Doppler effects (equations 30 and 28, respectively). The eventuall Doppler-hide results 
from the beaming pattern of the EC (star) emission, which - contrary to the discussed before 
SSC and EC(bl) processes - is maximised for small jet viewing angles (for a given [z/Ly] SJ/nifer 
one has [vS u ] ec ( s t ar ),br oc /+ T 2 ~ 5 2 ; equation 28). Comptonisation of the galactic infrared 
emission, as well as of the CMB radiation, is not expected to suffer such a decrease connected 
with the KN regime. For example, the Thomson regime condition for the EC(dust) process 
can be rewritten as v S yn,i4 < 60 B^^ST" 2 . For the standard nonrelativistic jet parameters, 
one can therefore put an upper limit [vS v ] ec (dust),br < 10~ 12 erg/s cm 2 at E ec (dust),br > 10 
GeV. At even lower photon energies the EC(CMB) emission dominates, proceeding well in 
the Thomson regime, with the maximum placed near £ ec {CMB),br ~ 0.1 GeV. The observed 
flux of this radiation is relatively low, [vS v ] ec (CME),br < 10~ 13 £>f 2 erg/s cm 2 . Detection of 
this emission from the FR I jets located further than a few Mpcs (i.e. further than Cen- 
taurus A) by future 7-ray telescopes is therefore possible only if the jet Doppler factor is 
significantly larger than unity, and/or ii B ^ B eq . On the other hand, large values of S are 
unlikely to occur in the FR I jets, although, as mentioned in section 2.1, the optical and 
X-ray jets can constitute the exceptions from this picture. 



3. Centaurus A and M 87 



The resulting inverse-Compton fluxes and break energies evaluated in the previous sec- 
tion depend on several free parameters. Because of relatively wide range of values covered 
by them, it is difficult to discuss 7-ray emission of the considered objects in general. Instead, 
a detailed analysis can be performed for each individual source. In this section, we discuss 
the large scale jet 7-ray radiative outputs of Centaurus A and M 87. One may note at this 
point, that both sources are quite unique. They are relatively nearby, what allowed in the 
past to study their jets at different frequencies and scales, with high spatial and spectral 
resolutions. The large scale jet in M 87 is the first one discovered at optical frequencies. 
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Both sources were among a few radio-loud AGNs known to posses the large scale X-ray jets 
before Chandra was launched. During the last decades, M 87 and Centaurus A were also 
being frequently observed at 7-rays, from sub-MeV to TeV photon energies, what resulted 
in positive detections or constraints on the flux upper limits. 



3.1. Centaurus A 



Centaurus A galaxy is a giant elliptical, which most probably mergered in the past with 
an other small and gas rich galaxy. With the distance D = 3.5 Mpc Centaurus A hosts the 
closest AGN, which is however difficult to be observed at IR and optical frequencies because 
of the obscuring dust lane. VLA studies (e.g., Burns et al. 1983; Clarke et al. 1992) show a 
double-lobed morphology of the Centaurus A radio source, with the jet and the counterjet 
of the FR I type and a kinetic power ~ 2 • 10 43 erg/s. VLBI observations (e.g., Jones et al. 
1996) followed the jet and the counterjet to the pc-scales, suggesting the jet viewing angle 
9 ~ 50° — 80°. Optical and IR observations of the Centaurus A nucleus (e.g., Bailey et 
al. 1986; Hawarden et al. 1993; Packham et al. 1996; Schreier et al. 1996, 1998; Marconi et 
al. 2000; Capetti et al. 2000) revealed the unresolved (< 1 pc) and variable central source 
with a small (~ 40 pc) nuclear disc and a possible extremely compact nuclear torus. No 
IR/optical emission connected directly with the jet non-thermal radiation was found (but 
see Joy et al. 1991; Schreier et al. 1996, 1998; Marconi et al. 2000). At the X-rays, the 
jet was observed in addition to the variable (at the time scales from minutes up to years) 
nuclear component (Schreier et al. 1979; Morini et al. 1989, respectively). Recently, Chandra 
telescope revealed details of the jet X-ray emission and presence of the X-ray counterjet 
(Kraft et al. 2002; Hardcastle et al. 2003), as well as of the whole Centaurus A system 
(Karovska et al. 2002). At higher energies, Centaurus A was observed by CGRO satellite, 
what allowed to detect the maximum of its nuclear 7-ray emission near ~ 0.1 MeV (Steinle 
et al. 1998). For e 7 > 0.1 MeV the 7-ray flux decreases, although at ~ 100 MeV photon 
energies it is still detectable (Sreekumar et al. 1999). CANGAROO observations of the 
extended (~ 14 kpc) region around Centaurus A nucleus put the upper limit on the emission 
at the VHE range, S(e 7 > 1.5 TeV) < 1.28 • 10' 11 ph cm' 2 s' 1 (Rowell et al. 1999). However, 
the discussed object was detected in the past at TeV energies, with the observed flux S(e 1 > 
0.3 TeV) ~ 4.4 • 10~ n phcm" 2 s" 1 (Grindlay et al. 1997), what made Centaurus A the very 
first (although not confirmed) detected extragalactic source of the VHE radiation. Recent 
upper limit suggests that this emission can be variable on a timescale of years, in analogy 
to the low- and high-states of activity known from the TeV blazar observations (but see also 
Harris et al. 1997, 2003, for the X-ray variability of the large scale jet emission in M 87). 
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Existence of a hidden BL Lac core in the Centaurus A nucleus, as expected in a frame- 
work of unification scheme, was widely discussed over the last decade (Bailey et al. 1986; 
Morganti et al. 1991; Hawarden et al. 1993; Packham et al. 1996; Steinle et al. 1998; Capetti 
et al. 2000). Most recently, Chiaberge et al. (2001) reconstructed broad-band spectrum of 
Centaurus A nucleus, from radio to 7-ray frequencies, and found spectacular similarities to 
the characteristic double-peaked blazar spectral energy distribution. The synchrotron com- 
ponent of this radiation was found to peak at far infra-red energy range, Vu,br ~ 10 12 — 10 13 
Hz, with the observed luminosity [vL v \ bl ,br ~ 10 41 — 10 42 erg/s. The inverse- Compton break 
energy was placed near ~ 0.1 MeV, with the observed power comparable to the synchrotron 
one. Chiaberge et al. fitted the SSC model to the multiwavelength Centaurus A nucleus 
emission, and found that all (except one) intrinsic parameters are similar to those of the 
low-luminous blazar sources. The only difference, as compared to the 'typical' BL Lac 
broad-band spectrum, was a small value of the required Doppler factor, Su ~ 1.6. Chiaberge 
et al. interpreted their results as the evidence for the jet radial velocity structure at pc-scales, 
consisting of a fast central spine surrounded by a slower boundary layer (see also Chiaberge 
et al. 2000). Assuming that the physical properties (i.e. electron energy distribution, mag- 
netic field intensity, etc) are the same in both jet components, the observed multiwavelength 
spectrum of the Centaurus A nucleus can then be regarded as a representation of a typi- 
cal low-luminous blazar emission, but originating within the slower jet boundary layer and 
therefore less beamed as compared to the 'classical' BL Lacs. It is consistent with the jet 
inclination ~ 70°. Most probably, Centaurus A observed at small angles to the jet axis would 
be therefore classified as LBL, with the observed luminosity ~ 10 45 — 10 46 erg/s and with 
the observed break frequency ~ 10 13 — 10 14 Hz. Note, that in such a case the Centaurus 
A nucleus is not expected to radiate at the VHE range (cf. Bai and Lee 2001). For the 
estimates below, we assume [uL v }u,42 ~ 0.3, i>bi,iA ~ 0.01, Tu ~ 10 and 5u ~ 1.6. 

In order to discuss the 7-ray emission of the large scale jet in Centaurus A, let us 
consider its brightest part in X-rays and at radio frequencies, the ensemble of the knots Al 
- A4. The X-ray observations (Kraft et al. 2002) suggest that for this region ~ 0.8, 
r\ ~ 0.4, £?_ 4 ~ 0.6, ax ~ 1-5 and L x ~ 4 • 10 39 erg/s. The synchrotron break frequency is 
unknown, as there is no observation of the jet synchrotron emission at IR/optical frequencies. 
However, the existing upper limits (Schreier et al. 1996) are consistent with the synchrotron 
spectrum resulting from the electron energy distribution given by equation 1, with ~ 0.5, 
Aa ~ 1 and the assumed v S yn,u ~ 1. Thus, below we put [vL^^n^,,. / Lx ~ (10 14 Hz/2.4 • 
lO^Hz) 1 -"*, implying for the Al - A4 region [vL v ] syTlj 42 ~ 0.2. Below we also take the 
galactic photons energy densities U star and Udust as estimated in section 2.2.3 7 . 



7 12 fim flux from Centaurus A is 23 Jy (Knapp et al. 1990). 
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For the considered parameters of Centaurus A, the IC fluxes and break energies in 
function of the jet viewing angle 9 are presented on Figure 1, for two different jet bulk Lorentz 
factors T = 2, 6. For small values of 9 the EC(star) and EC(dust) emissions dominate, while 
for large jet inclinations the SSC and the EC(bl) processes are the more important ones. In 
general, small jet viewing angles are excluded as they imply non-observed large VHE flux 
due to comptonisation of the starlight. However, for large jet inclinations the SSC radiation 
also seems to overproduce (comparing to the present upper limits) the VHE photons. For 
example, for the usually considered 9 = 70° and T < 2, the expected SSC energy flux is 
[vSvjss^br ~ 1CT 10 ergs -1 cm~ 2 at the observed break energy e SS c,br ~ 0.1 — 1 TeV. This 
suggest presence of the KN effects decreasing the SSC and the EC(star) emissions, what in 
turn indirectly suggests relativistic jet velocities (r > 2, Figure 2) and/or subequipartition 
magnetic field (B < 10~ 5 G, Figure 3). Alternatively, the unobserved synchrotron break 
frequency u syn ^ br can be higher than assumed here 10 14 Hz. 

For the moderate and large jet inclinations, the EC(bl) radiation is expected to be very 
strong, with the observed flux [vSy]^^),^- ~ 10~ n — 10~ 10 ergs -1 cm~ 2 and £ ec (u),br ~ 10 
GeV (for T = 1 — 10). This emission could be therefore easily observed in the future by 
the GLAST satellite. The EC(dust) emission is important only for 5 ^> 1, and otherwise is 
dominated by the EC(bl) radiation. Finally, comptonisation of the CMB radiation does not 
play a role in the 7-ray emission of Centaurus A unless 5 1 or B <C B eq (Figures 1 and 
3). 



3.2. M 87 

M 87 galaxy is a giant elliptical, placed in the center of the Virgo Cluster X-ray emission. 
With the distance D = 16 Mpc it hosts one of the closest AGNs, with spectacular and widely 
studied large scale structure. VLA observations of the M 87 radio source revealed one-sided, 
highly polarised and edge-brightened radio jet (e.g., Owen et al. 1989), in addition to the 
large scale cavities inflated in the intergalactic medium (Owen et al. 2000). A relatively 
high radio luminosity places the discussed object at the borderline between FR I and FR 
II radio sources, although the one-sided jet with a kinetic power ~ 5 • 10 44 erg/s displays 
FR I morphology. Detailed VLA and HST studies (Biretta et al. 1995, 1999, respectively) 
revealed apparent superluminal motions of the jet components at distances up to a few 
hundreds parsecs from the active center, while VLBI observations measured much lower 
velocities at parsec scales. This suggests a complex velocity structure of the M 87 jet, with 
highly relativistic bulk velocities over all of its length (r ~ a few) and a relatively small jet 
inclination to the line of sight (9 ~ 30°; e.g., Bicknell and Begelman 1996). IR observations 
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of the considered object (Whysong and Antonucci 2001; Corbin et al. 2002) give no evidence 
for the presence of an obscuring matter around the active center as expected for AGNs, 
consistently with a lack of the dusty nuclear torus or, eventually, with its small size (< 50 
pc). Optical observations of the M 87 nuclear region revealed a small disc (~ 100 pc) of 
ionised gas with a LINER emission line spectrum, fuelling a supermassive central black hole 
(e.g., Macchetto et al. 1997). The optical nuclear emission presents a smooth featureless 
continuum variable at time scales of months and connected with an unresolved (< 5 pc) 
central source (Tsvetanov et al. 1998). The large scale jet is very prominent and highly 
polarised in optical, what allows to perform its detailed radio-to-optical spectroscopy and 
polarimetry (e.g., Meisenheimer et al. 1996; Sparks et al. 1996; Perlman et al. 1999, 2001). 
At X-rays the jet and the nucleus were early observed by ROSAT satellite (e.g., Schreier 
et al. 1982). These observations allowed to construct broad-band radio-to- X-ray spectrum 
of the large scale jet emission (e.g., Biretta et al. 1991), to put constraints on the X-ray 
nuclear radiation and its variability (Harris et al. 1997) and, finally, to study interaction 
between the radio-emitting matter and the surrounding thermal gas within the extended 
M 87 halo (Bohringer et al. 1995, 2001). Recently, the X-ray jet radiation was studied in 
more detail by Chandra (Marshall et al. 2002; Wilson and Yang 2002; Harris et al. 2003). 
Besides several attempts in the past, till the last year M 87 was not detected in 7-rays, 
neither by CGRO satellite (Sreekumar et al. 1996) nor by ground-based detectors in the 
VHE range. Up to now, the existing upper limit on the TeV emission of the discussed 
object was S{e 1 > 0.21 TeV) < 1.2 • 10~ 10 ph cm" 2 s" 1 (Weekes et al. 1972). However, recent 
HEGRA observations resulted in positive detection of the VHE emission from the M 87 radio 
galaxy, with the observed flux S{e 1 > 0.73 TeV) ~ 0.96 ■ 10 -12 phcm~ 2 s~ x (Aharonian et al. 
2003). 

Similarly to the case of Centaurus A, presence of a hidden blazar in the nucleus of M 87 
was discussed previously by, e.g., Tsvetanov et al. (1998); Reynolds et al. (1999); Whysong 
and Antonucci (2001); Corbin et al. (2002). However, contrary to Centaurus A, a broad- 
band spectrum of the discussed central source in the IR-to-UV frequencies - and therefore 
parameters of the considered hidden blazar - is poorly constrained. One may note, that the 
radio-to- X-ray nuclear radiation is roughly similar to the large scale jet synchrotron emission, 
with a RO < 1 and a ox > 1 (Tsvetanov et al. 1998; Bohringer et al. 2001). This indicates, 
that the nuclear radio-to- X-ray continuum possibly results from the synchrotron emission 
of one electron population, as expected for the low-luminous blazars, with the observed 
synchrotron spectral break vu,br placed near optical frequencies. The observed synchrotron 
power is then ~ 3 • 10 42 erg/s (Tsvetanov et al. 1998). Furthermore, one may assume that, 
in analogy to Centaurus A, the observed M 87 blazar-like emission originates from a slower 
component of the blazar jet and, hence, that the Doppler factor of the discussed radiation 
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is close to unity. Here we take 1^,14 ~ 10 and [i/L u ]u, \i ~ 1 with the assumed Tu ~ 10 
and 5u ~ 2. The above values suggest that the considered M 87 blazar is situated on the 
borderline between LBL and HBL subclasses, and only large jet inclination preclude us to 
observe its 'classical' blazar properties (cf. Bai and Lee 2001). 

In order to evaluate the 7-ray flux of the M 87 large scale jet, let us consider its bright- 
est knot A, dominating radio-to-X-ray radiative jet output. Besides differences in the jet 
morphology as observed in different frequencies, multiwavelength studies (e.g., Biretta et al. 
1991; Meisenheimer et al. 1996; Perlman et al. 2001; Marshall et al. 2002; Wilson and Yang 
2002) allow one to construct broad-band spectrum of the considered region. They indicate 
R-i ~ 0.8 (at radio and optical; X-ray knots seem to be twice narrower), r\ ~ 1, £>_ 4 ~ 2.5 
and the synchrotron power ~ 10 42 erg/s. The exact value of the spectral break (and hence of 
the synchrotron break luminosity) depends on the detailed model of the synchrotron emis- 
sion assumed for the knot regions (Perlman et al. 2001; Marshall et al. 2002). However, 
the average optical spectral index «o ~ 0.9 (Perlman et al. 2001) indicates that v S yn,br is 
higher than 10 14 Hz. Below, for a rough estimate, we put v sy n, 14 ~ 10 and [vL v \ syn ^ ~ 0.3 
(Wilson and Yang 2002). We also take the galactic photons energy densities U sta r and Udust 
as estimated in section 2.2.3 s . 

A similar analysis as for the Centaurus A 7-ray emission suggests, that in the case of M 
87 large scale jet the EC(star) and EC(dust) processes should dominate production of the 
VHE photons (Figure 4). If the TeV flux detected recently by Aharonian et al. (2003) is in 
fact due to the knot A inverse-Compton emission, its most likely origin is comptonisation of 
the galactic stellar and circumstellar dust radiation. However, because of the expected KN 
effects, a detailed numerical analysis is required in order to find constraints on the jet bulk 
Lorentz factor and the magnetic field (cf. Figures 5 and 6). Possibilities that 5 > 1 and 
B < B eq could be additionally verified by the eventual detection of the EC(CMB) emission 
(Figures 4 and 6). The other IC processes for the M 87 jet are expected to be negligible 
because of the KN effects. 



4. Discussion and Conclusions 

In this paper we discussed the high energy 7-ray emission of the FR I large scale jets. 
We used the X-ray observations compiled with data from the other spectral bands in order 
to reconstruct energy distribution of ultrarelativistic electrons present in the considered 
objects. Our phenomenological approach to recover electron spectrum from observations 



12 /im flux from M 87 is 0.42 Jy (Knapp et al. 1990). 
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rather than from theories of particle acceleration is dictated by the fact, that such theories 
are still insufficiently developed to enable quantitative analysis. Next, we analyzed possible 
origin of the seed photons contributing to the inverse-Compton emission of the obtained 
electron spectrum, including nuclear jet radiation as well as ambient, stellar and circumstellar 
emission of the host galaxies, and discussed in detail spectral properties of the resulting 7-ray 
radiative output. The approach take into account effects connected with relativistic bulk 
motion of the emitting region, correcting earlier computations presented in the literature. 
We also clearly indicated limitation of the adopted approach due the Klein-Nishina regime 
effects. Other restrictions of the presented model are connected with only roughly evaluated 
galactic radiation fields and hardly known parameters of the hidden blazar radiation. Also, 
the effects of 7-ray absorption on CIB radiation are only shortly mentioned in this paper. 
However, even with all the aforementioned uncertainties, one can conclude that the detection 
of the high energy 7-ray emission by future ground-based and space telescopes from, at least, 
some of the FR I large scale jets is possible. Thus, the future observations will provide 
important constraints on the unknown jet parameters, like the magnetic field intensity and 
the jet Doppler factor, confirming or excluding possibilities of B <C B eq and 5^1 discussed 
in the literature. 

Our analysis provides constraints for 7-ray emission of the nearby FR I sources Cen- 
taurus A and M 87. For Centaurus A we predict measurable - by future 7-ray missions 
- fluxes at 10 GeV and 0.1 — 1 TeV photon energies due to comptonisation of the blazar 
radiation and the synchrotron self-Compton process, respectively. In the case of M 87 we 
show that recently detected VHE emission can result from comptonisation of the stellar and 
circumstellar infrared photons of the host galaxy. 

Possibility that the large scale jets - and not the active centers - in these two FR I 
radio galaxies can dominate production of high energy 7-rays was not discussed previously 
in the literature. Some other authors, e.g. Bai and Lee (2001), Protheroe et al. (2003) or 
Donea and Protheroe (2003), suggested and studied production of very high energy 7-rays 
in the nuclear regions of the considered sources by the small scale, blazar-like jets. This 
constitutes the main difference with our work. Unfortunatelly, the angular resolution of the 
I AC systems will not allow for separation of the kiloparsec scale jet 7-ray radiation from the 
eventuall core component. However, the core vs. the kpc-scale jet origin of 7-rays can be 
justified by studies of variability time scales. 
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A. The observed and intrinsic jet luminosities 

Let us consider an observer at rest in a frame K, and an emitting fluid comoving frame 
denoted by K'. Fluid velocity in the observer frame is (3c, the appropriate Lorentz factor is 
T, and the Doppler factor for a given inclination angle 6 = cos -1 fi with respect to the jet axis 
is 5 = [r (1 — (5 In K frame the radiation intensity is I u and the energy flux spectral 
density is S v = J I v dQ, where f2 = (cos -1 yU, 0). The total energy flux of the emission is 
S = JS V dv. 



As discussed by Lind and Blandford (1985) and Sikora et al. (1997), for a continuous 
(stationary) jet one has 

S» = ^5 J&W » ( A1 ) 

where j' u , is an intrinsic emissivity of the fluid (i.e. the one measured in K'), dV is the 
observed emitting volume element and D is the distance to the source. The observed 
isotropic luminosity is L = 4n D 2 S. Hence, with v — 8 1/ and dV = dV'/T, one obtains 

L = 4tt — — , (A2) 

where dL' /dfl' = J j' dV is the fluid intrinsic radiative power emitted in a given 
direction Q! = (cos -1 //', 0'), and j' = f j' v ,dv'. The total emitted power as measured 
in K' is then 

L'=ii^dn' = -L , (A3) 



on 1 5 3 

where the last equality holds for the intrinsically isotropic emission. 

For a moving single blob of emitting plasma (with the observed volume element trans- 
forming as dV = 5 dV) the observed energy flux spectral density is 

Su = ^ J & dV . (A4) 

Hence 

r)T' 

L = M 'mr • < A5 > 

and in the case of intrinsically isotropic emission one obtains 

/f)T' 1 
m dn ' = T> L ■ (A6) 
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Equations A3 and A6 specify the function g C j/ m b(T, 0) denned in section 2.2 for the 
two considered jet models. Note, that in both cases, the ratio between intrinsic power and 
comoving volume of the emitting region is (L'/V) = (L/V) 5~ 3 . Hence, the equipartition 
magnetic field is 

M*dWW*cS«*f ■ (A7) 

where v' eqmax = v e q,max/8 is the maximum synchrotron frequency considered in derivation 
of B eq . As a result, the equipartition magnetic field measured in the emitting region rest 
frame is related to the equipartition value computed for no beaming by relation 

B eq = B eq , &=l 5- b ' 7 . (A8) 

This equation corrects the formula A7 of Harris and Krawczynski (2002). 



B. External radiation fields 

Following the previous section, by K' we denote the inertial rest frame of the emitting 
matter of the kpc-scale jet. In the observer rest frame, K, the jet has bulk velocity (3c and 
the respective bulk Lorentz factor T. Let us consider radiation external to the jet. Energy 
density of this emission measured in K' is 

u' ext = - r [ iLtdV , (Bl) 

where I' is the jet comoving intensity of the external radiation and the integration is per- 
formed over the photons arrival directions. Defining next 

5 m = [r (1-/3 cos Vin)] -1 , (B2) 

where ip in is the angle between external photons and the jet axis, one has I' ext = I ext S~^, 
dQ! = dVt 5f n and thus 

ULt = \ [ h x t5- n 2 dn , (B3) 

c J AW 

• For emission isotropic in K frame in vicinity of the kpc-scale jet, one obtains 

u' ext = — ( C dn = u ext r 2 (i + ^) , (B4) 

C J47T 3 

where the last equality is due to U ext = Air I ext /c. 
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• In a case of the external emission illuminating the jet from its base one has ^ 
and hence 

1 f I ext dVt I ext AQ 1 



Uext ~cJ An , W~^~ (2fF ' ( } 

where I ext AQ is approximately the emission energy flux received by the stationary 
observer located at the jet axis at the distance r from the considered source of the 
emission. The isotropic luminosity computed by such observer is therefore L ext (0) = 
Air r 2 I ext AQ and 

T p _ L ext (0) 1 , 



C. Starlight emission 

Let us assume for the considered elliptical galaxies the spherical King-type stellar emis- 
sivity profile 

WC)« (i + C 2 ) _3/2 , (ci) 

where ( = r/r c , r G (0, r t ) is a distance from the galactic center, r c is a core radius and r t is 
a terminal boundary radius of the stellar distribution (see section 2.2.3). With the assumed 
emissivity one can find the appropriate radiative intensity, I s tar((, Ct), and next the mean 
intensity of the stellar radiation as a functions of r 

Jstar((, Ct) = J I star (C, Ct) dVt , (C2) 

where Ct = ft/fc- The analytical computations were performed by Tsai and Mathews 
(1995, see equation 26 therein). The energy density of the considered emission is simply 
U star (C, Ct) = (47r/c) J s tar(C) Ct) ■ Thus, knowing the central energy density of the starlight 
radiation, U sta r,c, one can evaluate U sta r at any distance from the galactic center as 

U s tar{Ci Ct) = U s tar, C 7^~Sr T\ ■ (C3) 

■Jstar{V, C,tJ 

For ellipticals with the bolometric luminosities ranging from Lboi = 1.14 • 10 44 erg/s to 
L bo i = 1.23 • 10 45 erg/s, Tsai and Mathews (1995) estimated U star:C = 1-2 • 10" 7 - 5.9 • 1(T 9 
erg/cm 3 , r c = 0.04—0.77 kpc and r t = 44.4—123 kpc, respectively (Table 1 therein). For these 
parameters, at the distance r = 1 kpc from the galactic center, the bolometric starlight energy 
density computed accordingly to the equation C6 ranges from 0.9- 10~ 9 erg/cm 3 up to 3- 10~ 9 
erg/cm 3 . Therefore, in this paper we take the average value £/ star (lkpc) ~ 10~ 9 erg/cm 3 . 
We also assume for simplicity an approximately isotropic distribution of the discussed stellar 
emission in the galaxy rest frame at the kiloparsec scale. 
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D. Anisotropy of the IC emission 

Let us consider the above defined rest frames K and K' . The emitting fluid radiates 
through synchrotron and inverse- Compt on processes, with intrinsic luminosities [v'L' v ,\ syn 
and [v'L' v ,] ic , respectively. The jet magnetic field energy density is U' B , and energy density of 
the seed photons as measured in K' is U' seed . The ratio of the observed isotropic luminosities 
due to both considered processes is 



r t 1 a (r a\ d i v 'K'^ d l u 'K> 

\vL v \ ic 4ng cj/mh (T,6)—Q$— 



[vLvUn A7rg cj/mb (r,e) 



(Dl) 



an' an' 



(Appendix A). Note, that in the above equation one compares the synchrotron and the 
inverse- Compton luminosities due to the electrons with the same energy 7 oc (v' / 1 B) 1 ^ 2 
and 7 oc (u' ic / 1 y' seed ) 1 ^ 2 (for example 7b r if break luminosities are considered). The intrinsic 
synchrotron power emitted in a given direction Q! = (cos -1 //, 0') by isotropic electrons with 
the energy distribution n' e {j) within an uniformly filled volume V is 

^^"KiW^n,,,^ , (D2) 

where \i\syn m( ? = § c °t U' b ^ 2 is a mean rate of electron energy losses due to the synchrotron 
emission. It is not difficult to show (see below), that in the Thomson regime considered in 
this paper an analogous expression representing comptonisation of anisotropic (in K') seed 
photons is 

— 0^} — = V [v j u , (S2 )\ ic =- — — — V |7(£2 ) \ ic mc , (D3) 

where |7(f2')|j C denotes electron cooling rate due to the inverse- Compton emission. For the 
seed photons antiparallel (+) or parallel (— ) in K' to the jet axis (if the seed photons are 
distributed isotropically around the relativistic jet, in its comoving frame their directions are 
almost antiparallel to the jet axis) one has 

l^')k ± ) = 2^7 2 (l±/i'f , (D4) 
while in the case of the isotropic (iso) distribution of the seed photons in the jet rest frame 

m%c (lso) = 4 - C ^U' seed 7 2 . (D5) 



Thus, the function f±/i SO (r, 6) defined in equation 16 is 

(iso) 



Mm -i it 1 *"? <t> . (D6, 
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Transforming it to the observer rest frame, with // = (// — — fifj,) one finally obtains 

/±/?so(r,^) as given by the equation 17 (see also Dermer et al. 1992; Dermer 1995). 

In order to find |7(fi')| ic for the cases discussed above, let us evaluate the appropriate 
inverse- Compton emissivity [j' u ,]i C - Here we write [v' i' v ,(£l')\ ic = mc 2 e' 2 h' ic (e', Q'), where the 
dimensionless photon energy e' = hu'/m e c 2 and the photon emissivity is 

n[ c (e',n>) = c J de' eed j dtt seed J j dQ' e 

(l-cosx')an' see ^'see d M' S ee d )n' e (^n' e ) . (D7) 

Here n' seed (e' seed , fl' seed ) is a seed photon number density, 774(7, fig) is the electron energy 
distribution, x' is the angle between the electron and the seed photon direction, a is the 
inverse- Compton scattering cross-section, and all quantities are measured in the emitting 
fluid rest frame K' (Dermer 1995). 



• For the seed photons antiparallel (+) or parallel (— ) in K' to the jet axis, assuming 
azimutal symmetry and characteristic seed photon energy (e' seed ), one can write 



n 'seed{ e 'seedi ^seed) ~ n seed $ ( 6 seed ( e 'seed)) $ (^seed ^ 1) ? 



(D8) 



where fi' 



seed 

rig(7, fi e ) = n' e (7)/47r and the Thomson scattering cross-section 

a = aT S[ e >- 7 2 e ' seed (1 - cos X ')} 6 [fi' - fi'J 
(e.g., Dermer 1995), one obtains 



(cos 1 fi' seed , (fi'seed)- Thus, with the isotropic electron energy distribution 



(D9) 



1 COt_ , 

2 4vr Useed 
din' e {rf)8 



i±n' 

€ ' ( e 'seed) 

7 - 



1/2 



<0(i±^') 



(DIO) 



This relation, with 7 = ^e'/{e' seed ) (1 5T/7) and U' seed = n' seed (e' seed ) mc 2 , leads to 
d[v' L' v ,] ic / dQ! given by expressions D3 and D4. 



For the seed photons isotropic in the jet comoving frame (iso), one can write 

n seed( e seedi ^seed) = n seed $ ( e seed ~ ( e seed)) 



(Dll) 
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Thus, with rig(7, Q' e ) = n' e (^)/ATi and the Thomson scattering cross-section 

4 



a = gt 5 



one obtains 



Kc(iso)( e 'i ^ ) 



^ 2 T ^seed 



1 C(J T , 

rn 

y d l n'Ja)& 



s [a - n'j 



3/4 

e ' ( e seed) 

7 - 



(D12) 



1/2 



4 <V' ) 

3 \^seedl 



(D13) 



This, with 7 = yV/f (e' seed ) and U' seed = n' seed (e' seed ) mc 2 leads to d[i/ L'M dW given 
by expressions D3 and D5. 



E. Break frequencies of the IC radiation 

Let us consider the break frequency of the large scale jet inverse-Compton emission. 
Accordingly to the formalism presented in Appendix D, in the jet comoving frame a seed 
photon with the characteristic (break) frequency v' seed is upscattered by the ultrarelativistic 
electron with Lorentz factor 7 to the frequency 7 2 v' seed (1 — cosx'). For the electron break 
Lorentz factor 7^ one obtains the observed break frequency of the inverse-Compton emission 

Vic(seed),br = Ibr Keed 5(1 ~ COS x') ■ (El) 

• For comptonisation of radiation distributed isotropically around the jet, one has v' seed = 
T v S eed and 1 — cosx' — 1 + ^ ■ Therefore 

Vec(seed), br = Ibr V seed T 5 (1 + fjf) = ^ V seed 5 2 || ^ ^ . (E2) 

• In the case of seed photons illuminating kpc-scale jet from behind v' seed = v see djV and 
1 — cos x' — \ — //. Hence 

Vec(seed), br = it "seed ^5(1-^')= 7br "seed & i 1 ~ lA { l + P) ■ ( E 3) 

• For comptonisation of the isotropic (in K') radiation one has v' seed = u seed /5 and the 
averaged W ic{seed) /v' seed ) = (4/3) 7 2 . Therefore 

4 2 

Vic(seed),br = ^Ibr^seed ■ (E4) 

Equations E2 - E4 specify the function h±/ iso (T, 9) defined in section 2.3. 
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Fig. 1. — High energy 7-ray emission of the Centaurus A large scale jet region Al - A4 
for different jet viewing angles 9, B = B eq and the bulk Lorentz factors T = 26. Other 
parameters are discussed in Section 3.1. The upper panels show the observed IC fluxes, 
while the lower panels indicate the respective observed photon break energies. The solid 
lines correspond to the EC(bl) emission, the dotted lines to the EC(dust) process, the short- 
dotted line to the EC(CMB) radiation, thick solid line to the EC(star) emission and, finally, 
the dashed and dash-dotted lines to the SSC radiation of a continuous jet and moving blob, 
respectively. The lines are presented only till the Klein-Nishina limit. 
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Fig. 2. — High energy 7-ray emission of the Centaurus A large scale jet region Al - A4 in 
motion of the bulk Lorentz factor T, for the jet viewing angle 9 = 70° and the magnetic field 
B = B eq . The other parameters and description of the presented curves are the same as in 
Figure 1. 
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Fig. 3. — High energy 7-ray emission of the Centaurus A large scale jet region Al - A4 in 
fuction of the magnetic field B, for the jet viewing angle = 70° and the bulk Lorentz factor 
T = 2. The other parameters and description of the presented curves are the same as in 
Figure 1. 
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Fig. 4. — High energy 7-ray emission of the M 87 knot A for different jet viewing angles 
9, B = B eq and the bulk Lorentz factors T = 26. The other parameters are discussed in 
Section 3.2. Description of the presented curves is the same as in Figure 1. 
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Fig. 5. — High energy 7-ray emission of the M 87 knot A in fuction of the bulk Lorentz factor 
r, for the jet viewing angle 9 = 30° and the magnetic field B = B eq . The other parameters 
and description of the presented curves are the same as in Figure 4. 
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Fig. 6. — High energy 7-ray emission of the M 87 knot A in fuction of the magnetic field B, 
for the jet viewing angle 9 = 30° and the bulk Lorentz factor T = 4. The other parameters 
and description of the presented curves are the same as in Figure 4. 



